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ABSTRACT 


The  purpose  of  this  work  is  to  analyze  and  evaluate  the 
effect  of  jamming  waveforms  on  both  coherent  and  noncoherent 
digital  communications  receivers.  Specifically,  random  processes 
are  utilized  as  jamming  models  in  which  it  is  assumed  that  the 
jamming  waveforms  have  been  produced  by  a  shaping  filter  driven 
by  white  Gaussian  noise.  Such  jamming  waveforms  are  then 
assumed  to  be  present  at  the  input  of  known  receiver 
structures  (in  addition  to-  the  signals  and  channel  noise 
normally  present) ,  and  optimum  jamming  waveform  spectra  are 
determined  for  different  receiver  schemes  and  modulation 
techniques . 

Graphical  results  based  on  numerical  analyses  are  presented 
in  order  to  demonstrate  the  effect  of  different  jamming 
strategies  on  receiver  performance.  In  order  to  quantify 
receiver  performance,  bit  error  probabilities  are  determined  for 
binary  modulation  systems  and  symbol  error  probabilities  are 
determined  for  M-ary  modulation  systems.  In  each  case,  the 
error  probabilities  are  functions  of  signal-to-noise  ratio 
(SNR)  and  j ammer- to-signal  ratio  (JSR) .  Results  show  that 
it  is  generally  possible  to  significantly  degrade  the 
performance  of  binary  as  well  as  M-ary  modulation  communica¬ 
tion  receivers  by  introducing  suitably  chosen  jamming  v’aveforms. 
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I. 


INTRODUCTION 


The  theory  of  statistical  signal  detection  and  estimation 
in  the  presence  of  additive  white  Gaussian  noise  is  widely 
described  in  many  textbooks  [Refs.  1,2,3].  Signal  detectors 
are  typically  designed  and  built  to  either  optimize  the 
receiver  output  signal  to  noise  ratio,  or  as  is  the  case  with 
digital  communications  receivers,  to  minimize  the  error 
probability . 

While  it  has  been  demonstrated  that  receivers  designed 
under  a  white  noise  interference  assiimption  tend  to  perform 
reasonably  well  even  when  the  interference  is  not  white 
[Ref.  4] ,  the  assumption  of  white  noise  interference  is  often 
invalid,  especially  when  the  receiver  must  operate  in  a  jamming 
environment . 

The  goal  of  this  thesis  is  to  analyze  the  vulnerability  of 
certain  digital  communications  receivers  designed  to  operate 
in  a  white  noise  interference  environment,  that  must  operate 
in  the  presence  of  jamming  also.  The  mathematical  model  of 
the  jamming  utilized  is  a  colored  Gaussian  noise  process 
whose  power  spectral  density  is  to  be  shaped  in  such  a  manner 
so  as  to  cause  a  large  increase  in  the  receiver  probability 
of  error.  While  it  is  not  always  possible  to  solve  certain 
spectral  shaping  optimization  problems,  it  is  possible  to 
postulate  techniques  that  intuitively  achieve  efficient 
utilization  cf  the  available  jammer  power. 


This  thesis  is  divided  up  as  follows.  In  Chapter  II,  we 
present  results  on  colored  noise  interference  effects  in 
coherent  M-ary  Phase  Shift  Keyed  (MPSK)  receivers,  and  receiver 
symbol  error  probability  in  the  presence  of  noise  and  jamming 
is  derived.  In  Chapter  III  we  analyze  and  determine  performance 
of  a  coherent  M-ary  Frequency  Shift  Keyed  (MFSK)  receiver 
operating  in  the  presence  of  noise  and  jamming.  Chapter  IV 
deals  with  non-coherent  Binary  Frequency  Shift  Keyed  (BFSK) 
signal  detection  in  the  presence  of  noise  and  jamming.  The 
performance  of  the  well-known  quadrature  receiver  is  analyzed 
under  dual  channel  and  single  channel  operation.  In  Chapter  V 
graphical  results  are  presented  and  discussed,  and  performance 
comparisons  are  carried  out.  The  conclusions  and  interpreta¬ 
tions  of  the  results  obtained  are  presented  in  Chapter  VI. 


II.  COLORED  NOISE  INTERFERENCE  EFFECTS  IN  COHERENT 
M-ARY  PHASE  SHIFT  KEYED  MODULATION 

A.  SIGNAL  DETECTION  IN  THE  PRESENCE  OF  COLORED  NOISE 

The  system  whose  performance  is  to  be  analyzed  is  described 
in  Fig.  2.1.  The  structure  shown  is  the  optimum  receiver  for 
recovery  of  MPSK  modulated  data',  in  the  presence  of  additive 
white  Gaussian  noise.  In  PSK  modulation,  the  source  (or 
modulator)  transmits  one  of  M  signals  s^(t)  ,  where 

i  =  1,2,...,N,  over  a  prescribed  time  interval.  Because  in 
transmissions  and  reception  these  signals  are  interfered  with 
by  noise,  at  the  receiver  one  observes  the  signal  r(t)  rather 
than  just  one  of  the  transmitted  signals.  Using  hypothesis 
testing  concepts,  we  say  that  under  hypotheses  ,  r(t) 

takes  on  the  form 


H.:  r(t)  =  VE  S.  (t)  +w{t)  +n  (t)  (2.1) 

1  1C 

0<t<T,  i=l,2,...,M 


where  for  M-ary  PSK  modulation 


(t) 


f2TTkt  2TT(i-l). 

V  Z/T  cos  ( - —  +  - - ) 

i  =  1 , 2  ,  .  .  .  ,  M 
k  is  an  integer 


,2.2) 


Here  W(t)  is  a  sample  function  of  a  white  Gaussian  noise 
process  of  Power  Spectral  Density  level  N  /2  and  n  (t)  is  a 


Scunple  function  of  a  colored  Gaussian  noise  process  having 


autocorrelation  function  K  (t )  .  We  assume  also  w(t)  and 

c 

n^(t)  are  statistically  independent  random  processes. 

The  receiver  of  Fig.  2.1  is,  as  previously  pointed  out,  a 
optimum  processor  (in  minimiam  error  probability  sense)  when 
n^(t)  =  0  .  The  analysis  that  follows  evaluates  the  effect 
of  on  the  performance  of  this  receiver.  Since 

may  represent  some  form  of  jamming,  the  error  probability 
expression  to  be  derived  can  be  used  to  determine  the  vulnera 
bility  of  such  a  receiver  to  colored  noise  jamming,  or 
conversely,  to  determine  the  colored  noise  spectrum  that  most 
effectively  causes  poor  or  inadequate  receiver  performance, 
namely,  high  error  probability. 

The  signals  S^(t)  ,  i  =  1,2,3,.. .,M  can  be  shown  to 
have  cross-correlation  coefficients 

T 

p..  =  /  s.(t)  S.(t)  dt  =  cos  (2.3) 

^  J  Q  J  ^  ^ 

i, j  =  1,2,  ...  ,M 

The  receiver  tal^es  advantage  of  the  fact  that  we  can  express 
the  S^(t)  functions,  i  =  1,2,...,M,  as  an  exact  (rather 
than  approximate)  expression  of  a  linear  combination  of  two 
functions  'p.  (t)  and  'p  (t)  .  In  other  words 


Cos  2^kt/T 

^t72 


(2.6) 


2  ( t ) 


Sin  2^kt/T 

kT/2 


where  k  is  an  integer. 

It  can  be  easily  shown  that 


---M—  '  "  ’  1-2 . « 


2-T(i-l) 

Sin  - r; -  .  1 

1*1  ' 


=  12  M 

9  ^  f  •  •  •  9  ^  ' 


%ji "  - 


S^(t) 


I  Sin  ^n^^^ 
n=i 


i  =  1/2,3,. ..,M 


(2.4) 


with 


Sin  "  ^  Si(t)  ii;n(t)dt-  n  =  1,2;  (2 

°  i  =  1,2,  .  .  .  ,M 

These  basis  functions  ii/i(t)  and  can  be  derived  via  a 

Gramm-Schmitt  orthonorir.alization  procedure  (or  almost  by 
inspection  in  this  case).  It  turns  out  that  ili(t)  and  V2(t) 
(which  must  be  orthogonal)  are  given  by 
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We  define 


9^  =  2TT(i-l)/M  i  =  1,2,.. .,M 


(2.10) 


and  assuming  equal  prior  probabilties ,  namely,  each  signal  is 
equally  likely  to  be  transmitted,  the  receiver  computes 


2 

I 


n=l 


i  =  1 , 2 ,  .  .  .  ,M 


(2. ii: 


and  makes  decisions  based  on  which  value  is  largest.  Thus 
with 


T 

r„  =  /  r(t)  ^^{t)  dt  n  =  1,2  (2.12) 

n  Q 


we  have 


a)j^(t)  dt 

r(t)  ^  ~ 


1,2,..,M  (2.13) 


and  using  simple  trigonometric  identities, 

=  V  Cos(e,  +-)  i  =  1,2,...,M  .  (2.14) 


Clearly 


V 


1/2 


(2.15) 


where 


T 

=  /  '  r(t)  ip,  (t)  dt  ,  (2.16) 

°  0  ^ 


T 

V  =  /  r(t)  i|/,(t)  dt 

®  0  ^ 


(2.17) 


and 


(2.18) 


B.  RECEIVER  PERFORMANCE 

Since  conditioned  on  any  hypothesis  i  =  1,2,...,M, 

V  and  V  are  Gaussian  random  variables,  we  can  obtain  the 
c  s 

statistics  of  the  appropriate  random  variables,  in  the  follow¬ 
ing  manner.  First,  we  have 


E;V  /H  .  ;■ 

c  : 


T 

E{  f  [)'£  S.(t)  +  w(t)  +  n  (t)]:^-,(t)  dt; 
o'  3  c  1 


_  T 

/E  /  S.(t)  V,  (t)  dt  =  .E  S.,  j  =  1,2,...,:-:  . 

3  -L  3 


(2.19) 


yn.  }  =  E{  /  [/E  s.  (t)  +  w(t)  +  n_(t)]ij;,  (t)  dt} 

»  J  n  -I  z 


=  /E  /  S  .  (t)  iijit)  dt  = 
n  J  ^ 


/E  S. 


(2.20) 


j  =  1,2,3,... ,M 


also 


T  2 

Var{V^/H.  }  =  E{[  /  [w(t)+n^(t)  ]i|/,  (t)dt]  } 

c  J  «  c-  J. 


E  {  /  /  [w  ( t )  +n_  ( t )  ]  [w  (  T )  +n_  ( T )  ]  ( t )  ip.  (  t  )  dtd  r  } 

0  0  c  i  ± 


N  T  T 

-^  +  /  /  K  (t-T)ip,  (t)ip,  (T)cit  dT 

0  0 


(2.21) 


Var(V^/H . } 
®  J 


T  2 

E{  [  I  [w  ( t) +n_  ( t)  ]  ip,  ( t)  dt]  } 
0  c  ^ 


T  T 


=  E(  /  /  [w  ( t) +n^  ( t)  ]  [w  (  T ) +n^  (t  )  ]  1^2  ( t)  ip^  ^  T )  dt  d~ 


0  0 


N  T  T 

+  I  J  K  ( t-r)  ip,  (t)  p  (  T  )  dt  dr 
^00 


(2.22) 


In  Appendix  A  we  demonstrate  that 


/  /  K  (t-T)  1/;,  (t)  I/;,  (t)  dt  dT 
0  0^  ^ 


*  A  A  ^ 

=  f  f  K  (t-T)  (t)  (T)dt  dr  =  a 

0  0  ^  ^  c 


(A. 7) 


SO  that  V  and  V  conditioned  on  H.  have  identical  variances, 
c  s  -  j 

Observe  also  that 


E{ [V^-E{V^/Hj }] [Vg-E{Vg/Hj }j/Hj } 


T  T 

E{^/  [w(t)+n^(t)  ]  ;i;^{t)dt  /  [w(T)+nj,  (t)  ]  :^2 


TT  N 


/f  -^-{(t-T)  'i;,  (t)  (T)dt  dl 

00  ^ 


+  /  /  K  (t-T) 'i/-|  (t)  (f)- (  t)  dt  dt 
0  0  ^ 


(2.23; 


VJe  can  observe  that  the  first  double  integral  in  Eq.  2.23 
is  zero,  so  that 


E-  [V^-E’ V  /H .  } ]  [V  -E{V  /H .  } 1/H  .  ; 
-  i.  c  ■  '  s  '  s  •  3  j 


it-:)  (t)  (  ■  )  dt  d' 


(2.24) 


We  demonstrate  in  Appendix  A  that  in  general  j,  9  is  not  zero 
so  that  V  and  V  conditioned  on  H-  may  not  be  uncorrelated. 


However  we  are  still  able  to  express  the  joint  probability 


density  function  of  V  and  V  by  using  the  general  form  [Ref. 

w  s 

5]  of  an  N-dimensional  Gaussian  random  vector  X,  namely, 


P^(X) 

where 


{27T 


N/2,,,1/2 


(2.25) 


m^  =  E{X} 


:2.26) 


and 


A  =  E^(X  -m^)  (X  -m^)^| 


(2.2’ 


In  our  case,  we  have  a  2-dimensional  problem  in  which  (see 
Eqs .  2.19  and  2.20)  , 


a 


,E 


3  =  1 , 2  ,  .  .  .  ,  M 


ana 


N 


—  +  ' 

2  'c 


2 

1.^ 


2 

1,2 


N 


2  c 


50  that 


(2.28: 


(2.29) 


12.30) 


It  is  simnie  to  snov; 


-ai_2/A 


(2.31) 


Thus 


Pv  ,V  /H,‘''c’V“j'  = 

C  S  j 


1  J  1  ^ 

— - exp  S  -2 

',-^,2 .  ^  V  -/ES.-, 

(2tt)-A  L  Ls  32  -J 


ta^j/A 


-“1,2/'' 


•“1,2/* 


(^+a^)/aJ  Lv^->'ESj2 


V  -  /ES  .  , 

C  :1 


:2. 32: 


with  j  =  1 , 2 , . . . ,M 


Now  we  need  to  obtain  from  this  probability  density  function  the 
joint  probability  density  function  of  V  and  n  conditioned  on  . 
This  type  of  transformation  [Ref.  6]  is  well  known  and  can  be 
used  here  to  obtain 


^V,n/H  ,  =  VP  /  'VCosr.VSinVH  ) 

^3  c  s  3 


+  VP,.  ,,  /„  (-V  Cos  n,-v  Sin  n/H  . )  ,  V  O, 

^c'V“3  ^  0  ~ 


(2.33) 


Using  the  probability  density  function  of  V^  and  V^  (Ecn.  2.3) 


yields 


%,VH  =  -= 

’  -J  ^2-)  i 


VCosn-vES.,  1  a  -b  VCosr-vES.,  '  I 

1  -  ^  ^ 
VSinn-^ES  .^J  [“b  a  jLvSin'  -vES]^  J  I 


[  ,  -VCosn-vES . ,  T  a  -b  -VCosr-vES.,  1 

^  -=exp  <^4  \,  V  ■  0,  (2.33) 

.2-:  -VSinn-vES._  -b  a  -VSinr— vRS._  f 

V  L  3.:J  L  JL  32JJ0  •  -. 
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a  =  (^  +  ol)/h 


(2.34) 


b  = 


(2.35) 


•  •  *  .*"1 


This  probability  density  function  can  be  expressed  in  the  form 


— -  exp  {“[a  (V  Cos  n  +  /ES^  ,  )  ^ 

'  2  2-1- 
(27r)  A 


+  a(V  Sinn  -/ESj2) 


•2b(V  Sin  n  -  '^'£8^2)  (V  Cos  n  -  /ES^^)]  } 


V  1  ""2  2 

-  exp  (--fa  (V  Cos  n  +  v'ES  . , )  -t-  a  (V  sin  n  +  >^ES  .  ^ ) 

- 5—  ^  3 1  ]  2 

(2^) ^A 


-  2b  (V  Sin  ■' :  +  >^'£5^2)  Cos  n  +  /ES^  j_)  ]  } 


(2.36' 


which  can  be  simplified  somewhat. 

Observe  that  the  exponential  of  the  first  term  simplifies  to 


a  [V^  E  -  AVk  E  (S  .  Cos  -  +  S  •  2  Sin  n)  ]  -  2b  [V^  Sin  n  Cos 


E(S.-,Sin’',  -(-S.^Cosn)  -i-ES.,S.- 

]1  j2  ]1  j2 


and  the  exponential  of  the  second  term  simplifies  to 


(2.37) 


a  [V  +  E  +  2V/E  (Sj^Cosn+Sj2Sinri)]  “  2b  [V  SinnCosn 

+  V/E(Sj;[_  Sin  n  +Sj2  Cos  n)  +  ££^^^3^2]  (2.38) 

We  can  now  group  certain  terms  together.  Observe  from  Eq.  2.8, 
Eq.  2.9  and  Eq.  2.10  that 


Cos  n  +Sj2  Sinn 


also 


Cos  Cos  n  -  sin  0^  Sin  n 
Cos  (0 j  +  n) 


S  j  Sin  n  +  S 


for  j  =  1,2,.. 


Cos  n  =  Cos  0j  Sin  n  ~  Sin  0 ^  Cos  n 
=  Sin  ( n  -  9 j ) 

.  ,M,  so  from  Eq.  2.36  we  have 


(2.39) 


(2.40) 


J 


'J  (2tt)  ^ 


^  exp{-|[a [V^+E-2V/ECos ( 9 . +"  )  ] 


2  - 

-  2b  [V  Sin  n  Cos  n  -  V>^  ESin  (n-  9^)  -E  Cos  0 ^  Sin  9  ^  } 


V 


^  (2i) 


12  — 

expi-y[a[V  +E+2V/ECOS (- .+n) ] 
^  D 


2  .  — 

-  2b[V  Sin’'Cos’’+  Vi'ESin(-';-9^)-E  Cos  9  ^  S  in  9  ^ 
for  V  >  0  and  0  r  ■ 


(2.41) 


Since 


sin (n  ~  9 j ) 


(2.43) 


=  -  Sin(r|-9j+'n’) 


we  have 


^  [exp{-i^[V^+E-2V/ECos  (0^+'^)  ] 


\/ (2tt) 


-  2b[V  Sin  n Cos  n-  V/ESin { n-0 j ) -ECos  e^Sin  9^]] 


.12—  2 
+  expi— 2[a[V  +E-2V/ECos  (  9  ^ +ri+'JT )  ]  -  2b  [V  Sin  n  Cos 


-  VvESin(n-9j  +  ^) -E  Cos  Sin9j]]}]  (2.44) 

for  V  ^  0  and  0  £  n 

It  is  apparent  from  the  range  of  n  that  the  two  exponential 
terms  can  be  replaced  by  a  single  term  with  v  ranging  from  0 
to  2  "  . 

Thus ,  we  have 


^  expi,-j[a  [V^+E-2VyECos  (-0  .+-  ) 


(2^^) 


-  2b[y  Sin'iCos"-VvESin(n-0j)-ECos0j  Sin-^j]  ;2.45) 


for  V  _  0  and  0  _  2".  The  probability  density  function  of 

conditioned  on  H_  is  obtained  via  integration  of  P,^  iV,’  ,  H^)  , 

J  r  '  J 

name Iv  . 


.2.46) 


Returning  to  our  decision  rule,  (Eq.  2.14)  ,  recall  that  we 
decide  baaed  on  which 


=  VCps(e^+n)  i  =  1,2,. ..,M  (2.14) 


is  largest. 


So,  if  is  the  true  hypothesis,  then  a  correct  decision 


is  made  if 


V  Cos  (  3  j+n  )  >  V  Cos  (3^+n)  ;  i  =  1 , 2  ,  .  .  .  ,M 

i  j 


(2.47) 


Since  Cos  x  is  maximum  when  ;x;  is  minimum,  we  see  that  if 
Hj  is  the  true  hypothesis,  a  correct  decision  is  made  if 


^  =  ^'2 . M 

i  ^  j 

Now  from  Eq .  2.9  we  know  that, 

=  27(j-l)/M 

So  Eq .  2.54  is  satisfied  for  ■"  in  the  region 


(2.48) 


]  M 


+  — 
j 


(2.49) 


Thus,  the  probability  of  making  a  correct  decision,  given  that 
H  is  the  true  hvoothesis,  Pr-c.H.  ,  is  aiven  bv 


Pr{c/Hj } 


(2.50) 


-fl  +i 


/ 


If  we  make  the  variable  change 


S  =  n  +  0- 


(2.51) 


Then  Eg.  2.50  becomes 


tt/M 


Pr{c/H.}  =  /  P  /H  (B-e./H  )dS 

J  -tt/M  j 


(2.52) 


Now  from  Eq.  2.45  and  Eq.  2.46  we  have 


P,/„  (n/H,)  =  / 


-  exp{-4la  1V^+E-2V/E  Cos 


2b  [V^  Sin  n  Cos  n-V/ESin ( n-e j ) -E  Cos  9^  Sin  8^ ] j  }dV  (2.53) 


0  <  n  <  2" 


so  that 


/  —  ^  -  -  ■  exp{-j[a  [V^+E-2Vv  E  Cos 


°  V^(2^)^A 


-  b [V^Sin2 (3-9 . ) -2V/ESin{3-2- •  ) -E  Sin 2-j ] ]  ^  dV 


2.54) 


. 


and  £q.  2.52  now  becomes 


TT /M  “  _ 

Pr{c/H. }  =  f  f  - - -  exp{-4[a[V^+E-2V/ECos  6 

-  b[V^Sin  2(B-0j)-2V/ESin(6-2ej)-E  Sin 26j] ]  }  dV  dB  (2.55) 

Since  the  hypotheses  have  been  assiuned  to  be  equally  likely, 
we  have 


Pr  (c) 


M 

M  ^  Pr{c/Hj} 
M  j=l  ^ 


(2.56) 


so  that 

P  =  1  -  Pr'.c; 

e 

=  1  -  -  [  I  I  — ■■■  -  -  exp{-i[a[V  +E 

^  j  =  l  -tt/M  0  ^^2.)  2, 

—  2  — 

-  2V.'E  Cos  2-b  [V  Sin  2  (3-0^  ) -2V/ESin(i-29j  ) -E  Sin  2-5  j  ]  ]  :-dV  di 

[2. SI) 

Observe  that  if  colored  noise  is  not  present,  then  from  Equations 
2.24,  2.30  and  2.35,  1  =  (N^/2)^  and  b  =  0,  so  that  Eq.  2.57 
simplifies  to  the  well-known  expression  for  the  performance  of 
the  M-PSK  receiver  operating  in  the  presence  of  additive  Vv’GN. 


That  is. 


-  2v/e  Cos  B+E]  ]  }dV  dB 


(2.58) 


where  in  Eq.  2. SI,  the  dependence  on  the  index  j  disappears 
when  b  =  0.  While  Eq.  2.57  yields  a  mathematical  result  on 
the  performance  on  the  M-PSK  receiver  in  the  presence  of 
WGN  and  colored  noise  jamming,  its  further  analysis  represents 
a  separate  project  in  itself.  Not  only  must  Eq.  2.57  be 
optimized  for  energy  constrained  jamming  but  also  it  must  be 
evaluated  when  the  jamming  spectrum  takes  on  some  simple  forms. 
For  this  reason,  no  effort  has  been  made  to  further  develop 


the  above  results. 


III.  COLORED  NOISE  INTERFERENCE  IN  COHERENT  M-ARY 
FREQUENCY  SHIFT  KEYED  MODULATED  SYSTEMS 

A.  SIGNAL  DETECTION  IN  THE  PRESENCE  OF  COLORED  NOISE 

The  structure  of  the  demodulator  whose  performance  is 
to  be  analyzed  is  shown  in  Fig:-  3.1.  This  receiver  is 
known  to  be  optimiam  for  deciding  with  minimum  probability  of 
error,  which  one  of  M  different  signals  forming  an  orthogonal, 
equal  energy  set  received  in  additive  white  Gaussian  noise  was 
actually  transmitted.  The  problem  analyzed  here,  can  be 
stated  as  follows:  A  waveform  r(t),  received  in  the  interval 
(0,T),  contains  one  of  the  M  signals,  Sj^(t),  i  =  1,2,...,M, 
with  equal  probability,  as  well  as  white  Gaussian  noise  w(t) 
of  Power  Spectral  Density  level  N^/Z  and  colored  Gaussian  noise 
n^{t)  having  autocorrelation  function  K^(t).  The  signals  are 
orthogonal  with  energy  e.  That  is 


1] 


S^{t)Sj(t)dt  = 


1  =  ] 
i  j 


(3.1) 


The  decision  rule  used  by  the  receiver,  is  to  choose  S^(t) 


as  the  transmitted  signal  if  G^  is  a  maximum,  where 


G,  =  f  r(t)  S,  (t)  dt  i  =  1,2,. ..,M 


(3.2) 


Fiyure  3.1.  WSK  im.-oivor 


While  this  is  an  optimum  test  (in  minimum  probability  of  error 


sense)  in  the  absence  of  the  colored  noise  n  (t) ,  the  analy- 
sis  of  the  next  section  is  carried  out  in  order  to  determine 
the  effect  of  n^(t)  on  the  receiver  performance.  Since 
will  typically  be  inserted  in  the  channel  by  an  unfriendly 
jaunmer,  it  is  reasonable  to  assume  that  n(t)  and  n^(t)  are 
statistically  independent  random  processes. 

B.  RECEIVER  PERFORMANCE 

Since  is  the  output  of  the  ith  correlator,  and,  condi¬ 
tioned  on  any  hypothesis,  is  a  Gaussian  random  variable, 
we  can  obtain  the  appropriate  conditional  statistics  that  allow 
determination  of  P^,  namely  the  receiver  error  probability. 

Thus 

T 

E-.GVH  ;  =  E{  /  [S.  (t)  +  w{t)+n,.(t)  ]S  .  (t)dt;- 

]  i  0  c  j 

T 

=  J  S.  (t)S.  (t)dt  =  ^  (3.3) 

0  ^ 


and 


T  2 

VarGVH.;  =  E-[  |  [  w(  t )  +  n  ( t )  ]  S  _  ( t )  dt  ] 

J  1  g  c  j 


31 


.  ' 


TT  N 

GVH.j  =  //  <S(t-T)  +  K  (t-T)  ]S.  (t)S  .  (T)dtdT 


NT  T  T 

=  I  S .  (t)S  .  {t)dt  +  //  K  (t-T)S.  (t)S.  (T)dtd 

^  Q  J  J  OQ  ^ 

N  TT 

=  ^  e  +  II  K  {■t-T)S.  (t)S  .  (T)dt  di  (3. 

2  c  -  3  3 


Define 


so  that 


Var{G ./h  } 
J  ^ 


N 

o 

2 


C, 


+ 


2 

0 

c,3 


(3. 


Observe  furthermore  that 


E-  [G.-E{G^/H^}]  [G^-E-'Gj^/H^}]/H^  } 


E'  /  [w(t)^n^{t)  JS.  (t)dt  '  [w(r)+n^(’:)  ]S,  (:)d^ 

0  ^  J  0  c  K 


T  T  M 


T  T 


j  ^  -  (t-:)S  .  (t)S.  (-)dtdr  j  K  (t-  :  )  S  .  ( t)  S,  (  ' )  dtd- 


E{  [Gj-E{Gj/H^}]  [Gj^-E{Gj^/H^}]/H^} 

NT  TT 

=  /  S.(t)S,  {t)dt  +  //  K  (t-T)S  .  (t)S,  (T)dt  dT 

^0^  OqC  3  < 

NTT 

=  -2  e  6.J.  If  K  (t-T)S.  (t)S.  (T)dt  dT  (3.7) 

^  00  ^  3  . 

As  can  be  seen  from  Eq.  3.7,  due  to  the  presence  of  the  colored 
noise,  the  random  variables  {Gj/H^^}  are  not  uncorrelated. 
However  we  will  show  that  for  MFSK  with  signal  frequencies 
that  are  sufficiently  separated,  the  integral 

TT 

II  K  (t-T) S  .  {t)S,  (T)dt  dT  (3.8) 

00'^  3  k 

vanishes  for  j  k,  so  that  the  random  variables  are  indeed 
uncorrelated . 

Thus,  conditioned  on  the  G^  are  statistically  indepen¬ 

dent.  Assume  now  that  Sj^(t)  is  transmitted  and  =  x.  Then 
the  conditional  probability  of  a  correct  decision,  Pr-c/Hj^, 


G-  =X'  becomes 


P  ic/H^  t  —  X  }  =  P  {G^  <  X  /  •  «  .  ,  ^i-1  ^  t  •  *  •  t  <  x/H^  ,  G 


n  P{G,  <x/H.  ,G.  =x} 
k=l  ^  ^  ^ 


Introducing  a  change  of  variable. 


k/^i 

Now,  since 


we  have  that 


p{c/H. }  =  /  P{c/K.  ,G.=x}  - = - exp{-(x-e)  /2  (^+o  .)}dx 

-®  N  ^  c,i 


so  that  using  Eq.  3.11,  we  obtain 


X/ 

00  M 

pic/H^}  =  /  n  / 

-oo  k  =  l  -00 

kj^i 


exp  j- (x-  e)^/2(-^£+a^^^  )}dx  (3.12) 


Assume  now  for  convenience  that  M  is  odd,  and  express  the  M-ary 
FSK  signals  as 


S^(t)  =  A  Cos  (uj^  +  (i- (M+l)/2)  Auj)  t  0  <  t  <  T  (3.13; 

i  =  1 , 2  ,  .  .  .  ,  M 


so  that 


T  -j.^t  * 

f  S.  (t)e  dt  =  /  S.  (t)p(t)e 

0  _oo 


=  F{S[  (t)p(t)  :•  i  =  1,2,. ..,M  (3.14) 


where 


p(t)  = 


1  0  <  t  <  T 


P(aj)  = 


u)T/2 


0  otherwise 


and  S|(t)  is  just  S^(t)  with  -oo  _<  t  £  ».  Thus 


F{S!_(t)p(t)  } 


00 

=  ^  f  SI  (w-A)P{A)dA 

2  T  ■'  1 


where  with 


(M+1),, 
(1  -  — 2 - ' 


we  have 


S[{'a))  =  tA  [  5  {-j- )  +  5  (uj+ (to^+’xi^)  )  ] 

for  i  “  i^2/***fM  • 


Thus 


.%  w'^  •  V'* 


(3.15) 


(3.16) 


(3.17) 


(3.18) 


F-  S|  (t)  p(t)  '' 


-]..T/2  ^ 


f  tA  [1  (oj-v-  )  +  •.  (  1  "Fs  - 


a'. 


-j(--'^  -..-.)T/2  Sin(.-..  -...)T/2 
AT  Cl  Cl 

=  -^le  — ^ — 


--..)T/2 
c  1 


+  e 


—  j  (  -hxi  ■ )  T/2  sin  (  ju+'j-.  "Hji./ . )  T/2 

C  X  Cl 


.+•.,  +’..)T/2 
c  1 


(3.19) 


For  convenience,  let 


■  I  ""1  ■■?  "’J  "T'^T- 


-.»  V 


L  (co)  =  e 


So  that  Eq.  3.79  becomes 


-jajT/2  Sin  a!T/2 
a)T/2 


(3.20) 


F{S!_(t)p(t)}  =  F{S^(t)}  =  S^(aj)  = 


AT 

=  -:r-[L(oj-aj  ~Cjj ■  )  +  L(arhja  +aj- )  ]  (3.21) 

2  Cl  Cl 

i  =  1 , 2  ,  .  .  .  ,  M 


Let  us  examine  now  the  correlation  coefficient  ^  ij  ' 


11 


/  S.(t)S-(t)dt  =  /  ACos  (oj  +tjj  •  )' tACos  (  ^  ,  )  t  dt 

,  1  J  n  c  j.  Cl 


A^T 


Sin  (u  . -0) .  )  T  Sin  (2w_+(jo  . +w  . )  T 


1  1 
(cu^-ojj  )  T 


c  1  : 


(2w  +a) .  +(jJ  .  )  T 
i  1 


(3.22; 


If  we  assuiae  that  -ai^T  >>  i,  then  the  second  term  in  Eq.  3.2  2 
vanishes  and  we  have 


11 


0 


S^  (t) Sj  (t) dt 


A  T  Sin(i-j)  1a.T 
2  (i-j)A,..,T 


3.23; 


In  order  to  have  orthogonal  signals  we  need  at  least  1,.T  = 
or  equivalently  1^  =  -/T.  Normally,  we  will  have 


k-/T 


(3.24) 


where  k  is  a  large  integer,  so  that  Pj_j  ~  0  i  7^  j  •  Thus, 

from  Eq.  3.23  and  Eq.  3.24, 


)..  =  /  S.  (t)  S.  (t)dt  = 

-‘■J  A  J 


e  =  A  T/2  i  =  j 


i  7^  j 


;3. 25: 


2 

From  Eq.  3.5  it  appears  however  that  the  term  a 

c,  1 

independent  of  i.  Nevertheless  Eq.  3.12  becomes 


LI  ~ 

0^2 
■^+0  , 
2  c  ,k 


P-.c/H^}  = 


k=l  -o 


V  2  c ,  1 


expi -(x-s)  ^2(^+0^  . 

2  c ,  1 


(3.26) 


o  ^  2 
2  c ,  1 


(3.27) 


Then 


(3.28) 


So  that  Eq .  3.26  becomes 


p{c/H^}  =  /  n  / 

— oo  Jc=l  — OO 


- e  '  dz 

/27T 


1  -n  /2  , 

X  - e  ^  dn 


Finally , 


Pic}  = 


1  ^ 

k  I  P':c/Hi} 
i=l 


or  equivalently 


/n 


M  »  M 


P'-C;  =  ^ 


V  -  1  ■■ 

1=1  k=l 


erfc. 


/_  _ 

C\/f--c,K 


-.2/2 


dn 


We  must  focus  on  the  J  ,  factor.  Observe  from  Eq.  3." 

c ,  k 


E-  [G^-E-G./hV  •]  [G,  -E;G,  /H^  }]/H,  ; 

J  J  -i-  ^  K.  X  X 


T  T 

.  +  ■'  K  ( t- '  )  S  .  ( t )  S,  (  : )  dt  d  ■: 

00  ^ 


(3.29) 


(3.30) 


(3.31) 


that 


and  the  second  term  becomes 


If  K  (t-T)S.(t)S.  (T)dtdT  =  If  K  (t-T)S’(t)S’(T)dtdT 
0  0  J  ^  J  < 


2tt 


J 


S  (ui)  S (-co)  S'  (co)  du) 

C  J  K 


(3.32) 


It  has  been  shown  in  Appendix  B  that  Sl(-cj)  and  S/(aj)  are 

3  ^ 


essentially  freque.'icy  disjoint,  therefore  Eq.  3.32  is  zero 
for  j  ^  k.  For  j  =  k,  we  have  (using  Eq.  3.21) 


G,k 


^  (uj)  lS^(^.o)  l^do) 

— oo 


=  “2^  ‘  ^  1 L  (ui-oj^-ojj^)  +  L(arHo^+cu^)  j dijj  (3.33) 


If  we  define 


T  f  2. 

Ij^  =  2?r  ^  |L(u)-w  -Wj^)  +  L(aj+io^+'j0j^^)  ]  dw 

>*oo 

k  =  1,2, ... ,M 


then,  with  e  =  A  T/2  (Eq.  3.25) 


2 

'c ,  k 


e  Ij^  k  =  L  ,  2  ,  .  .  .  ,H 


(3.34) 


(3.35) 


Thus  from  Eq.  3.31 


(3.40) 


S  (w)  =  2'n’K  'l  [6((jo+(jj  +0).  )  +  6  (u)-co  “oi .  )  ] 

c  i=l  c  i  Cl 


Thus,  the  colored  noise  consists  of  equally  weighted  "tones" 
at  the  signal  frequencies.  Therefore,  Eq.  3.34  becomes 


Ir  =  ^  /  2tK  J  [6  +  6  ]  |L(lij-'o^~loj^)+L(uj+co^+cjj^)  ;  ^du. 

-00  i=l 


M 


,2. 


M 


=  TK  ^  [  lL(-2(jj^-ojj^-aj^)  +L(cjj^-aj^)  1^  +  |L(a)^-L0j^)  +L(2tj^+r^+Li)^)  1^]  (3.41) 


i=l 


Since  oj^  is  typically  large,  we  can  justify  the  statement  that 


the  terms  involving  2oj^  are  negligible  small,  so  that. 


M 


S  TK[  1^  +  iL(aj^-(^)  1^] 


M 


2TK  )  i  L(uj.-oj^)  j  =  2TK 
i=l'  ^ 


M  /Sin(u.-co^)T/2\ 


M  /Sin  'i-k)AwT/2  \  2 
2TK  [ 


\  (i-k)/l,.T/2 


) 


3.42: 


With  l'...T/2  =  my  where  m  is  large,  we  have 


I,  =  2TK  for  i  =  k 
k 


(3.43) 


Ue  can  impose  a  constraint  that 


2?  J  S^(a))du) 


K  J  [6  (ij+ui|,+cOj^)  +  6  ((Jj-oj^-ojj^)  Idw 

-oo  i=i 


rhen 


(3.44) 


Furthermore 


K  =  P  ./2M 


(3.451 


=  2T 


(3.46) 


I,  TP  . 

k  _  n  j 

£  Me 


(3.43) 


and  since 


TP^j  =  jammer  energy  and 
e  =  signal  energy, 


this  implies  that  TP  ./e  =  JSR.  We  have  therefore  that 

n] 

Eq.  3.39  becomes 


^e  ^'^+-|[o.54.jsr- SNR/M  )] 


1  -n  /2  , 

—  e  dr 


=  1-  j 


i  [l-erfc.  (  n*  )]  ^e" 


-'/2 


(3.47) 


Observe  that  for  JSR  =  0,  Eq.  3.47  is  identical  to  the  well 
known  formula  for  the  perfoirmance  of  the  receiver  of  Fig.  3 
under  MFSK  modulation. 


IV.  NON-COHERENT  BINARY  FREQUENCY  SHIFT  KEYED 

SIGNAL  DETECTION  IN  THE  PRESENCE  OF  COLORED  NOISE 

A.  THE  QUADRATURE  RECEIVER,  EQUIVALENT  FORMS  AND  RECEIVER 

PERFORMANCE  IN  THE  PRESENCE  OF  WHITE  GAUSSIAN  NOISE 

In  this  section,  a  short  presentation  of  the  basic  princi¬ 
ples  of  statistical  communication  theory  that  lead  to  the  design 
of  the  well-known  quadrature  receiver  is  undertaken.  Basic 
results  that  are  useful  in  the  sequel  are  presented  only,  since 
the  details  have  been  worked  out  in  numerous  textbooks  (see 
[Ref.  7]  for  example). 

Consider  a  binary  digital  communication  system  model  in 
which  one  of  two  signals,  SQ(t)  or  S^(t),  with  energy  E^  and 
E^,  respectively,  is  received  in  the  time  interval  (0,T).  At 
the  receiver,  white  Gaussian  noise  with  zero  mean  and  spectral 
density  Nq/2  is  added  to  the  signal.  The  actual  received 
signal  r(t)  takes  on  one  of  the  two  forms,  namely 

r(t)  =  ^E~  S^{t)  +  n(t)  0  t  £  T,  i  =  0,1  (4.1) 

The  likelihood  ratio  test  which  operates  on  r{t)  in  order  to 
choose  which  one  of  the  two  hypotheses  is  believed  to  be 
the  true  one,  namely 

r(t)  =  yE~  S^{t)  +  n(t),  0  _  t  _  T,  (4.2) 

1  =  0,1 


45 


I 


i-  •  u 


t 


I 


■ 

I 


is 


exp{- i  /  [r(t)  -  1^  Sj_(t)  ]  )dt 

A(r(t))  =  - - j—  J  7  (4.3) 

exp{~  /  [r(t)  -  ^  S^(t)  ]  }dt 
o  0 

where  ;  is  a  threshold  whose  value  depends  on  the  decision 
criteria  used.  This  test  can  be  applied  to  any  communication 
problem  involving  transmission  of  known  signals  S^(t)  and  S2(t). 
One  such  example  is  the  well-known  BFSK  modulation  scheme. 

Ine  nrcblem  of  interest,  which  is  a  slight  modification  of 
BFSK  modulation  problems,  involves  signals 


vE^  S^tt;  =  A  Sin(uj^t+?^)  i  =  0,1 

0  <  t  <  T 


(4.4) 


where  the  phases  i  =  0,1  are  statistically  independent 

random  variables,  uniformly  distributed  over  the  interval 
(0,2~),  and  the  amplitudes  A  are  known  and  equal.  It  turns 
out  that  the  test  specified  by  Eq.  4.3  can  be  modified  to 
account  for  the  random  phases  by  using  conditional  proba¬ 
bility  densities. 

The  details  of  the  procedure  have  been  worked  out  in  Reference 
8.  It  can  be  shown  that  when  the  signals  are  given  by  Eq.  4.4, 
the  test  of  Eq.  4.3  becomes 


■  (r (t) ) 


I^(2Aqi/N^) 


(4.5) 


46 


where 


2 


■  T  -|2  T  2 

/  r(t)Sin  ojj^tdt  /  r(t)Cos  1  ,  k=  0,1 


(4.6) 


and  Iq(’)  is  the  modified  Bessel  function,  defined  by 


Io(x)  = 


I 

n=0 


2n 


2^"(nJ)^ 


■ 


2tt 


^xCos ( 6+a) 
e  da 


(4.7) 


For  minimum  error  probability  decision  criterion,  the  decision 
rule  of  Eq.  4.5  assuming  equal  prior  probability  of  trans¬ 
mitting  SQ(t)  or  S^(t),  is  to  choose  if 

I^(2Aq3_/N^)  >  l^(2AqQ/N^)  (4.8) 

or  equivalently,  to  choose  if 

2  2 
^1  -  % 

Otherwise  is  chosen.  (Observe  that  I  (x)  is  a  monotonicallv 
o  o 

increasing  function.) 

The  receiver  structure  that  implements  the  test  of  Eq. 

4.8  is  shown  in  Fig.  4.1.  Another  (equivalent)  form  of  the 
receiver  of  Fig.  4.1  is  shown  in  Fig.  4.2,  involving  a  combina¬ 
tion  of  matched  filters  and  the  envelope  detectors.  The 


AV' 


receiver  of  Fig.  4.2  is  completely  equivalent  to  the  receiver 
of  Fig.  4.1. 

The  evaluation  of  the  performance  of  the  receiver  has  been 
worked  out  in  Reference  9  and  is  given  by 

,  -E/2N 

Pg  =  ^  e  °  (4.9) 

2 

where  E  =  A  T/2  is  the  average  signal  energy.  If  we  now 
define  the  signal  to  noise  ratio,  (SNR)  as 

SNR  =  E/Nq 

we  obtain  the  simple  result 

Pg  =  J  exp{-  SNR/2}  (4.10) 

B.  RECEIVER  PERFORMANCE  IN  THE  PRESENCE  OF  COLORED  NOISE 

The  receiver  presented  in  Section  A  is  optimum  in  minimum 
probability  of  error  sense  when  operating  in  a  white  Gaussian 
noise  interference  environment.  In  this  section  we  analyze  the 
vulnerability  (probability  of  error)  of  the  quadrature  receiver 
in  the  presence  of  an  additional  additive  npise  that  is  modeled 
as  colored  and  Gaussian,  having  autocorrelation  function 
( T  )  .  (We  denote  i^^,(t)  as  this  additional  colored  noise). 

The  problem  can  then  be  restated  as  follows.  Under 
hypotheses  H. ,  i  =  0,1,  r(t)  takes  on  the  form 


r(t) 


Sj_{t)  +  w(t)  +  n^(t) 


i  =  0,1 
0  <  t  <  T 


(4.11) 


where 


S^(t)  =  /E  S[(t)  =  A  Sin(w^t  +  ())^)  i  =  0,1  (4.12) 


In  order  to  determine  the  effect  of  n  (t)  on  the  receiver 

c 

probability  of  error,  we  evaluate  the  statistics  of  the  random 
2 

variables  k  =  0,1,  where,  as  defined  by  Eq.  4.6, 


T 

2 

T 

=  /  r(t)Sin  0).  tdt 

_  0  . 

+ 

/  r(t)Cos  to.  tdt 

-  0  . 

k  =  0,1 


Thus,  conditioned  on  H^,  i  =  0,1 


.  (t)  +  w  (t)  +  n^  (t)  ]  Sin  w,  t  dt 

1.  C  K 


+  /  [S.  (t)  +w(t)  +n  (t)]Sin  co,  t  dt 

n  1  C  k 


2  2 

^i,k  ^  ^i,k 


i  =  0,l  k=.  0,1 


(4.13) 


Observe  first  that  the  integral 


/  S  ■  (t)  sin  0),  t  dt  =  /  A  Sin  (cj.  t  +  <j) .  )  Sin  w,  t  dt 

0  ^  0  1  1  K 


AT 

2 


Sin(w.-w,  )T/2 

(aj^-Wj^)T/2  Cos  [  (w^-a)j^)T/2 


sin  ((jj . +tjj,  )  T/2 

(w.+w,  )T/2  Cos[  (aj^+^^^)T/2  +  ^ 

X  iC 


i  =  0,1 
k  =  0,1 


(4.14) 


If  we  now  assume  that 


(ojt-(jj_)T  =  2mTr  and  (a)-,+a)  )T  =  2l7r 
i  o  -  1  o 


(4.15) 


we  have  that  Sin(w^  +coj^)T/2  =  0,  for  i  =  0,1,  k  =  0,1.-  Thus 


T 

/  S-  (t)  Sin  0),  tdt 

Q  1  K 
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Sin(Wj^-cjj^)T/2 

(Wi-oj,^)T/2 


Cos  [  1/2  +  ^>j_]  0^ 


i  =  0,1  k  =  0,1 


(4.16) 


where 


1  if  i  =  k 
0  if  i  k 


(4.17) 


By  arguments  similar  to  the  above. 


T 

j  S.  (t)Cos  '^tdt 


AT 

9 


Sin((^^-:.^)T/2 
' -  )  r/7 


> 


Sin[  (■...-^)T/2  + 


(4.18) 


are  Gaussian 


Now  conditioned  on  41 .  ,  i  =  0,1,  the  X.  ,  and  Y.  , 

1.  ^  ^  ^  f  ^ 

random  variables,  so  it  is  possible  to  obtain  the  conditional 
probability  density  function  of  qj^,  k  =  0,1.  Thus 

E{X.  ,/h.  d).  }  =  /  S.(t)Sin  oj.  t  dt  i  =  0,1  k  =  0,1  (4,19) 


and 


r  T 


Var-.X^^^H^,  =  E 


/  (w(t)-tn^(t)  )Sin  -iij^tdt 


rp  ip 


j  /  EK(t)  w(t)  +  n  (t)  n  (t)  }Sin  lax  tSin  ...Tdt  di 
00  c  c  K  K 


T  N 


1  y 

j  ~  Sin  u)  tdt  +  /]  K  (t-T)Sin  ^t  Sin  ^  idt  dT 

2  iC  r\f\  ^  XX 


tt 


00 


’w  ^'c,k  ^  ^ 


(4.20) 


where,  assuming  that  2xj^T  >>  1 


_2 
'  w 


and 


N  T 

~  j  Sin^.;,tdt 
2  Q  X 


N  T 
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k  =  0,1 


(4.21) 


T  T 


K  (t-:)Sin  t  sin  .j,  :dt  d-  k  =  0,1 


:4 . 22: 


Similarly 


•V-j.'vt' 


E  {  Y.  ,0  }  =  /  S.  (t)  Cosuj.tdt  i  =  0,1  k  =  0,1  (4.23) 

1/ivll  Q  X  K 


T  N  TT 

Var(Y^  k/^i''^i^  ~  ^  Cos*^ ojj^tdt  +  //  Kjt-T)Cos  Wj^t  Cos  ‘^"dt  d* 


A  2^2 
=  a  +  o  , 
w  c,k 


k  =  0,1  i  =  0,1 


since  it  can  be  demonstrated  that 


T  T  T  T 

I  /  K  (t-T)Sin  cjutSin'Aj,  dtdx  =  //  K  (t-T)Cos  co,  t  Cos  <-o, 'dtdi 

O'Q  c  K  Od  k  k 

k  =  0,1 


Finally 


E  [x.  -EiX,  ,/ri.,o.r  Y.  ,-E{y.  ,/H. 

jL  i,k  i,k/  1  1  JL  i,k  i,]c  1,  1  J  1 


T  T 

Ej  /  [w(t)  +  r^(t)  iSin  -^tdt  j  [w( t)  +  r^( t)  ] Cos  d- 


^Sin  t  Cos  -vtdt  +  i  '  K  (t-r)Sin  t  Cos  'dt  d- 
0  2  k  K  00  c  ^  k 


i  =  0,1 


(4.24) 


(4.25: 


(4.26) 


It  can  be  shown  that  these  two  integrals  are  zero  so  that 


and  ^  are  conditionally  uncorrelated,  and  therefore 
independent  since  they  are  Gaussian  random  variables.  Now 
define 
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A 


2 


k  =  0,1 


and 


2 

^k 


o 


2 

c  ,k 


k  =  0,1 


(4.27) 


(4.28) 


so  that  the  conditional  density  functions  for  k 


=  0,1  are, 


1  (  (q’+A’  )' 


(^> 


where  u(*)  is  the  unit  step  function,  and 


=  E^{X  /H  ,:■}+  I Y  /H  ; 

OO  0,0  o  o  o,o'  o  o 


Using  Eq.  4.16  and  Eq.  4.18,  we  have 


,  ,AT  ^  ,2  ,AT  ^  .  ,2 

oo  =  'o^  ^  (^Sin 


(^)2 
'  2 


Also 


?(qyHi,:i)  =  — y  exp 


4.29: 


(4.30) 


(4.31, 


where 


+  E^{Y^^Q/H^,(j)^}  =  0  (4.32) 


due  to  the  result  of  Eqs.  4.16,  4.17  and  4.18  for  i  7^  k. 
Therefore 


P(qyHi,(})il 


exp{-  qy2a^}u(q^) 


(4.33) 


Furthermore 


(  ^^^^0,1^  .  1 

'V^i^o,i  \ 

I  2al  °  ' 

(4.34) 


where  again  due  to  Eqs.  4.16,  4.17  and  4.18, 


=  0  (4.35) 


so  that 


- 2  exp^  -qV2j  'u(q,  ) 

2a 

^1 


(4.36) 


Finally 


P(q’/H  , :  )  = 


(4.38) 


We  now  have  the  statistical  information  needed  to  compute  the 
probability  of  receiver  error  P^.  Assuming  that  each  hypothe¬ 
sis  has  equal  prior  probability,  we  have 

^e  =  I  P^^l  7^- ^"^1  ‘"^o 

Observe  that 


pCqi-^o 


/  P(qi>q„/H^.q„  =  p)p  (p/H^><Jo 

oo 


(4.40) 


where 


OO 

P^^l  "  V^o'%  =  /  P(q3_/H^)dq^ 

0 

I 
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(4.41) 


Since  the  conditional  probabilities  functions  are  not 
dependent  on  the  individual  phases,  that  is 


P(q^/H^)  =  P(q^/H^,0^)  i  =  0,1 


(4.42) 


and 


P(qyHJ 


P(q 
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0,1 


(4.43) 


we  can  rewrite  the  conditional  probability  functions  in 
the  following  form 
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Similarly 
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so  that 
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Using  now  Eq.  4.44  and  Eq.  4.46  we  have 
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For  convenience,  let  =  AT/2  and ,  recalling  u(o)  =  1,  c 
Eq .  4.52  becomes 
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(4.53) 


(4.54) 


(4.55) 


(4.56) 


(4.57) 


Now  the  integral  itself  yields  1,  since  it  is  the  integral 
of  a  probability  density  function.  Therefore 


O  '  20_ 


T 


(4.58) 


Similarly  for  Eq.  4.51 


p{q.-<i  <o/ii.}  =  /  r  / 

—00  L 
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(4.59) 


Observing  first,  the  quantity  in  braclcets  can  be  expressed  as 
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Letting 


=  AT/20]_ 


(4.61) 


and  making  a  change  of  variable 


we  have  that  Eq.  4.60  becomes 


/  2  2,., 

“  -(x  +a,)2 

1-  /  xe  I^{aTX)u(x)dx 


p/a. 


O  1 


1  -  Q(a^,p/a^) 


(4 


where  Q(*,-)  is  the  well-known  Marcum  Q  function  [Ref.  10]. 
Therefore  Eq.  4.59  becomes 
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From  Reference  11,  the  integral  of  Eq.  4.64  becomes 
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Thus  Eq.  4.59  becomes 
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Now  using  Eq.  4.58  and  Eq.  4.59  in  Eq.  4.39,  we  have  that 
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Recalling  that  e  =  AT/2,  using  Eq.  4.55,  we  have 
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From  Eq.  4.67 
must  maximize 
noise  power. 


it  is  clear  that  in  order  to  minimize  P  ,  we 
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where 
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As  an  example,  consider  the  case  where  the  power  spectral 
density  of  the  jammer  is 


=  ttP 
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Under  this  assiamption  it  has  been  shown  in  Appendix  C,  that 
Eq.  4.69  becomes 
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where  oj-  and  P  are  the  frequency  and  the  power  of  the  jamming 
waveform,  respectively.  It  has  also  been  demonstrated  in 
Appendix  C  that  Eq.  4.7  0  is  maximum  at  cu, 
that 
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and  Eq .  4.68  now  becomes 
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Thus  Eq.  4.67  becomes 


2  ,  2 

where  SNR  =  (A  T/2)/N  and  JSR  =  P  T/(A  T/2),  represent  signal 

o  ^ 

to  noise  ratio  and  jamming  to  signal  ratio  respectively. 

Observe  that  with  JSR  =  0,  Eq.4.73  becomes  identical  to 
Eq.  4.10.  This  result  is  appealing  because  for  the  case  of 
no  jamming,  the  receiver  performance  should  be  identical  to  that 
of  a  receiver  operating  in  white  Gaussian  noise  interference 


C.  RECEIVER  PERFORMANCE  IN  THE  PRESENCE  OF  WHITE  GAUSSIAN 

NOISE  UNDER  SINGLE  CHANNEL  OPERATION 

In  Section  B,  we  have  analyzed  the  performance  of  the 
quadrature  receiver  in  the  presence  of  white  and  colored 
Gaussian  noise.  Results  were  specifically  obtained  when  the 
colored  noise  interference  was  a  single  frequency  jammer.  Sup¬ 
pose  now  that  the  quadrature  r-eceiver  experiences  a  single 
frequency  interference  which  corresponds  to  one  of  the  signal 
frequencies,  say  Since  the  receiver  makes  binary  decisions 

based  on  whether  q^^  >  q^  or  vice  versa,  the  presence  of  the 
interference  at  frequency  will  cause  q^  to  be  greater  than 
q^  most  of  the  time  creating  decision  errors  nearly  50%  of  the 
time. 

In  order  to  prevent  this  type  of  situation  from  arisini^ , 
the  receiver  can  turn  off  the  affected  channel,  or  equivalently, 
make  decisions  based  only  on  the  output  of  the  other  channel, 
that  is,  based  only  on  the  size  of  q^^ .  In  this  section  the 
performance  of  the  quadrature  receiver  is  analyzed  assuming 
white  Gaussian  noise  only  interference,  and  that  decisions  based 
on  only  one  channel  output  are  made. 

Assuming  that  the  receiver  bases  decisions  only  on  the 
size  of  q^^,  the  decision  rule  now  becomes 


Recall  from  Eq .  4.6  that. 
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The  probability  of  error  is 
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and  assuming  that  PiH^}  =  P{H^}  =  1/2  then  Eq.  4.76  becomes 
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and  in  Section  B  we  found  that 
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where  e  =  AT/ 2  and 
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Thus,  from  Eqs .  4.19,  4.46,  and  4.47,  we  obtain 
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Observe  however  that  a  threshold  of  y  must  now  be  defined. 
Clearly,  a  threshold  that  minimizes  should  be  chosen.  This 
can  be  done  by  solving  dP^/dv  =  0. 

Thus 
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so  that  solution  of  dP^/d :  =  0  yields  an  implicit  solution  for 
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Suppose  now  that  Yq  the  solution  of  Eq.  4.80  for  a  given 
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Letting  y  =  q^/a^,  Eq.  4.81  becomes 
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so  that  defining 
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we  have 
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so  that  the  threshold  setting  equation  (Eq.  4.80)  becomes 
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and  Eq.  4.82  simplifies  to 
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The  receiver  performance  indicated  by  Eq.  4.87  is  compared  to 
that  of  an  incoherent  BFSK  receiver  that  utilizes  both  channels 
for  its  decisions.  (See.  Eq.  4.10.)  The  result  of  this 
comparison  is  presented  in  Chapter  V. 

D.  RECEIVER  PERFORMANCE  IN  THE  PRESENCE  OF  COLORED  NOISE 
UNDER  SINGLE  CHANNEL  OPERATION 

In  this  section,  we  analyze  the  performance  of  the  quadra¬ 
ture  receiver  under  the  assumption  of  single  channel  operation, 
as  described  in  the  previous  section.  Here  however,  it  is 
additionally  assumed  that  a  jamming  signal  is  present,  whose 
energy  is  concentrated  around  the  frequency  (Observe  that 


69 


the  channel  whose  output  is  has  a  passband  around  Thus 

a  jammer  concentrating  its  energy  around  would  significantly 
affect  the  output  q^.  Consequently,  turning  off  or  ignoring 
q^  would  make  sense  under  these  circumstances.  Hence,  the 
single  channel  operation  being  considered  here. ) 

Our  decision  rule  continues  to  be 
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Observe  that  due  to  the  presence  of  a  jammer 
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As  shown  in  Appendix  C, 
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when  the  jaunmer  is  concentrated  at  frequency 
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so  that  the  probability  of  error  is 
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we  have 
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Observe  that  with  JSR  =  0,  Eq.  4.92  becomes  identical  to  Eq. 
4.87,  as  must  be  the  case. 

Furthermore  if  the  frequency  separation  is  such 

that  T/2  >>  1  or  T/2  =  mir ,  where  m  is  an  integer 

then,  SSQ  becomes  very  small  or  zero  so  that  the  effect 
of  the  presence  of  the  jamming  is  negligible.  The  numerical 
results  obtained  from  Eq.  4.92  are  very  similar  to  those 
obtained  from  Eq.  4.87  as  demonstrated  in  greater  detail  in 
Chapter  V. 

Recall  that  the  threshold  is  obtained  from  the  solution  of 
Eq.  4.86,  namely 


(4.92) 


I^(  /^^(2SNR) )  =  e^^^  (4.86) 

However  if  our  goal  is  to  set  a  threshold  that  minimizes  , 
for  the  case  being  considered  here,  we  can  solve  for  an  opti¬ 
mum  threshold  setting  by  minimizing  Eq.  4.92  with  respect  to 
r_„.  If  this  procedure  is  carried  out,  we  obtain  the  threshold 
setting  equation 
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(4.93) 


While  this  result  is  intuitively  appealing,  a  practical 
problem  arises  in  that  in  most  cases,  the  receiver  does  not 
know  the  operating  JSR  value,  hence  a  threshold  could  not  be 
set. 

Fortunately,  computer  evaluations  carried  out  using  both 
Eq.  4.86  and  Eq.  4.93  to  set  the  threshold  have  demonstrated 
that  the  resulting  with  thresholds  set  by  Eqs.  4.86  and 
4.93  are  almost  (and  for  all  practical  purposes)  identical. 


GRAPHICAL  RESULTS 


A-  GRAPHICAL  RESULTS  FOR  COLORED  NOISE  INTERFERENCE  IN 

COHERENT  M-ARY  FREQUENCY  SHIFT  KEYED  MODULATED  SYSTEMS 

In  Chapter  III,  the  performance  of  the  MFSK  receiver  in 
the  presence  of  white  and  colared  noise  was  derived.  This 
mathematical  result  is  used  now  to  evaluate  and  graphically 
display  receiver  performance  under  the  presence  of  white 
noise  only  and  under  the  presence  of  white  and  colored  noise 
interference . 

Results  are  presented  sequentially  for  values  of  M  =  2 ,  4, 

8,  and  16  on  the  performance  of  the  M-ary  FSK  receiver  for 
white  noise  as  the  only  source  of  interference  as  well  as  for 
various  conditions  of  colored  noise  powers  in  addition  to  the 
normally  present  WGN  interference.  The  performance  results 
for  the  M-ary  FSK  receiver  presented  in  this  section  in  terms 
of  the  probability  of  error  are  shown  as  the  SNR  changes , for 
specified  values  of  JSR.  Some  representative  results  are 
summarized  in  Tables  5.1,  5.2,  5.3  and  5.4.  Figures  5.1  through 
5.4  include  the  performance  of  the  M-ary  FSK  receiver  when  the 
transmitted  signal  is  interfered  by  white  noise  only,  namely, 

JSR  =  0.  This  makes  it  possible  to  evaluate  the  effect  of  the 
jamming  on  the  receiver  in  comparison  to  the  case  in  which 
WGN  is  the  only  source  of  interference.  These  results  have 
been  obtained  by  evaluating  Eq.  3.47. 
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B.  GRAPHICAL  RESULTS  FOR  NON-COHERENT  BINARY  FREQUENCY 
SHIFT  KEYED  SIGNAL  DETECTION  IN  THE  PRESENCE  OF 
COLORED  NOISE 

In  Chapter  IV,  the  performance  of  the  quadrature  receiver 
operating  in  the  presence  of  white  and  colored  noise  was 
derived.  The  mathematical  results  are  now  used  to  evaluate 
and  graphically  display  receiver  performance  under  various 
conditions  of  signal  and  noise  powers. 

First,  results  are  presented  for  the  case  in  which  white 
noise  is  the  only  source  of  interference.  This  yields  the 
well-known  probability  of  error  curves  for  the  standard  quadra¬ 
ture  receiver  for  non-coherent  BFSK .  These  are  presented  in 
Fig.  5.5,  along  with  a  corresponding  plot  of  the  probability  of 
error  of  the  quadrtature  receiver  in  which  only  one  channel 
output  is  used  to  make  binary  decisions. 

Additionally,  the  performance  of  the  quadrature  receiver 
operating  in  the  presence  of  white  and  colored  noise  is  evalu¬ 
ated  under  dual  channel  and  single  channel  operation.  Under 
single  channel  operation,  it  is  assumed  that  the  colored  noise 
jamming  concentrates  its  energy  around  one  of  the  FSK  operating 
frequencies,  and  that  the  receiver  is  able  to  make  a  determinis¬ 
tic  as  to  which  "channel  is  being  jammed"  so  that  the  outputs 
of  this  channel  are  ignored  in  the  process  of  making  decisions. 
Evaluations  are  carried  out  using  receiver  thresholds  that  are 
dependent  as  well  as  independent  of  jamming  power  levels. 

(Both  cases  are  considered  separately.)  The  performance  of 
the  quadrature  receiver  in  the  presence  of  noise  and  the  jamming 


waveform  described  in  this  section  in  terms  of  the  probability 
of  error  is  calculated  as  the  SNR  changes  for  specified  values 
of  JSR.  Some  important  results  are  summarized  in  Table  5.5 
for  JSR  =  0  and  in  Tables  5.6-5.10  as  JSR  takes  on  values 
of  0.0  db,  5.0  db,  10.0  db,  15.0  and  20.0  db,  respectively. 

In  Figure  5.5  the  performance  of  the  standard  quadrature  receiver 
and  the  single  channel  operation  of  the  quadrature  receiver  is 
plotted  when  the  transmitted  signal  is  interfered  by  white 
noise  only.  The  theoretical  performance  of  the  standard 
quadrature  receiver  is  calculated  from  Equation  4.10,  and  the 
performance  of  the  quadrature  receiver  under  single  channel 
operation  is  calculated  from  Equation  4.87. 

In  Figures  5.6-5.10,  the  performance  of  the  standard 
quadrature  receiver  and  the  quadrature  receiver  under  single 
channel  operation  with  the  threshold  dependent  as  well  as 
independent  of  the  jamming  power  level  is  plotted  when  the 
transmitted  signal  is  interfered  by  white  noise  and  by  the 
jamming  waveform  having  Power  Spectral  Desnity  given  by  Equation 
C.7.  Each  of  the  figures  corresponds  to  a  specific  value  of 
JSR  as  shown  in  the  headings.  The  performance  of  the  standard 
quadrature  receiver  is  calculated  from  Equation  4.73.  The 
theoretical  results  for  the  single  channel  operation  of  the 
quadrature  receiver  with  a  threshold  that  is  independent  of  the 
jamming  power  (Eq.  4.86)  is  calculated  from  Equation  4.82,  and 
Equation  4.92  is  used  to  compute  performance  of  the  same  re¬ 
ceiver  when  the  threshold  is  dependent  on  the  jamming  power 
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level.  As  pointed  out  in  Section  D  of  Chapter  IV,  the  proba¬ 
bility  of  error  calculated  from  Equation  4.92  with  the  thres¬ 
hold  set  by  Equations  4.86  and  4.93  show  almost  identical 
results . 

Tables  5.5  through  5.10  demonstrate  that  the  performance 

of  the  quadrature  receiver  under  single  channel  operation  is 

unaffected  by  changing  values  of  JSR.  This  is  due  to  the  fact 

that  for  (jj,  and  T  values  used  in  the  simulation,  the  value 

o  1 

of  SSQ  term  in  Eq.  4.92  is  identical  to  zero.  Thus  in  order 
to  demonstrate  the  effect  of  the  jammer  on  the  receiver  under 
single  channel  operation,  the  value  of  the  jamming  frequency 
has  been  allowed  to  vary  from  all  the  way  up  to 
Thus,  in  place  of  the  SSQ  term  as  defined  in  Eq.  4.92,  we  use 
the  modified  term 


SSQ  = 


Sin  {to  T/2 

(Wj-coj^)  T/2 


to  <  tu  .  <  to, 

0-3-1 


The  results  of  these  modifications  are  presented  in  Fig.  5.11 
and  Fig.  5.12  where  the  probability  of  error  of  the  receiver 
given  by  Eq.  4.93  is  evaluated  for  JSR  =  5  db  and  JSR  =  10  db , 
respectively,  where  the  jamming  frequency  (jo^)  is  allowed  to 
take  on  values  o-  =  (which  corresponds  to  the  results  given 
by  Eqs.  4.91  and  4.92  without  modification),  and  values  of 
=  3  (  +  -q)/4  and  Some  of  the  important  results 

obtained  are  summarized  in  Tables  5.11  and  5.12  for  JSR  = 


5  db 
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12.  Performance  of  the  quadrature  receiver  single  channel 
operation  for  different  jamming  frequencies  ana 
JSR  =  10  db 
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TABLE  5.11 

PERFORMANCE  OF  THE  QUADRATURE  RECEIVER  SINGLE  CHANNEL  OPERATION 
FOR  DIFFERENT  JAMMING  FREQUENCIES  AND  JSR  =  5  DB 
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TABLE  5.12 

PERFORMANCE  OF  THE  QUADRATURE  RECEIVER  SINCLE  CHANNEL  OPER 
FOR  DIFFERENT  JAMMING  FREQUENCIES  AND  JSR  =  10  DD 
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VI.  CONCLUSIONS 

The  analysis  carried  out  in  this  thesis  presents  the 
application  of  concepts  derived  in  statistical  communication 
theory,  specifically  in  the  theory  of  signal  detection  under 
the  assumption  of  colored  noise  interference.  The  performance 
of  digital  receivers  in  terms  of  probability  of  error  is 
determined  when  the  receivers  operate  in  the  presence  of  white 
and  colored  Gaussian  noise.  Three  techniques  are  examined 
separately,  one  for  MPSK  modulation,  another  for  coherent 
MFSK  modulation  and  the  last  one  for  (incoherent)  BFSK 
modulation . 

The  mathematical  model  of  the  jamming  waveform  proposed, 
consists  of  colored  Gaussian  noise  of  different  spectral  shapes 
and  power  content. 

For  MPSK  modulation,  a  mathematical  result  on  the  performance 
of  the  (coherent)  receiver  in  the  presence  of  WGN  and  colored 
noise  jamming  was  derived.  The  complexity  of  the  result  along 
with  the  many  possible  trade-offs  involving  spectral  shapes, 
power  levels  and  frequencies  of  operation  made  it  impossible 
to  address  in  this  thesis  the  issue  of  optimum  jamming  strate¬ 
gies  for  MPSK. 

For  MFSK  modulation  results  on  the  effect  of  the  coherent 
receiver,  were  derived.  A  simple  assumption  was  made  on  the 
spectrum  of-  the  jamming.  By  assuming  that  each  signal  frequency 


t  < 


<  ( 


was  interfered  with  a  tone  subject  to  a  total  jaitiming  power 

constraint,  the  receiver  was  evaluated  for  different  values 

of  SNR,  JSR,  and  M.  The  results  demonstrate  that  this 

form  of  jamming  can  be  quite  effective  or  that  significant 

increases  on  can  be  achieved  even  at  low  JSR  values, 
e 

For  the  case  of  BFSK  modulation,  the  quadrature  receiver 
was  analyzed  under  two  conditions  of  operation,  standard 
operation  and  single  channel  operation,  in  the  presence  of 
colored  noise  jamming  with  different  power  levels.  The  single 
channel  operation  was  introduced  as  a  method  for  mitigating  the 
effect  of  a  single  tone  jammer  at  one  of  the  carrier  frequen¬ 
cies.  When  no  jamming  is  present,  single  channel  operation 
performs  slightly  worse  than  standard  receiver  (both  channels) 
operation.  However,  in  the  presence  of  jamming,  single  channel 
operation  is  superior  to  standard  operation  because  the  receiver 
is  capable  of  eliminating  much  of  the  jammer  energy  and  its 
effect  by  ignoring  the  output  of  the  jammed  channel  during 
single  channel  operation.  As  pointed  out  in  Chapter  IV,  the 
effect  of  the  jamming  waveform  on  the  receiver  under  single 
channel  operation  depends  strongly  on  the  jamming  frequency 
chosen.  For  the  single  channel  operation,  it  was  assumed  that 
the  jamming  is  present  at  one  of  the  two  signal  frequencies,  and 
that  the  receiver  turns  off  the  channel  affected.  Thus,  deci¬ 
sions  are  made  based  only  on  the  output  of  the  unaffected 
channel.  However,  if  under  this  condition  of  operation  the 
jamming  changes  its  frequency  in  such  a  way  as  to  "move 


closer"  to  the  frequency  of  the  unaffected  channel,  it  has 
been  demonstrated  that  the  receiver  probability  of  error  in 
creases  as  Wj  approaches  the  frequency  of  the  unaffected 


channel . 
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APPENDIX  A 

DETAILED  INVESTIGATION  OF  THE  VARIANCES  OF  Vp  AND  Vq 
CONDITIONED  ON  HYPOTHESES  Hj 


Let 


T 

I 


''  ip^(t)  0  <  t  <  T  Z  =  1,2 

•< 

0  otherwise 


(A.l) 


and 


F{,i;^(t)  }  I  =  1,2 


(A.  2) 


Thus 
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I  ~  If  ip'^(r)dt  dT 

'  — <30-X 


OO  00 


/  /  (co)e^‘^^^~^^dwip'y(t)ii^J(T)dt  dT 


*00  *oo 


/  .■)  \  w/  V  «  ' 
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^  /  S  (a))d(j 


(A.  3) 


where  K  { r ) 


S_(oj  ) 


i  J./  ^ 


e"^  (co-2TTr/T)T/2 


(a)-^)T/2 


+  g-j((i^2TTn/T)T/2 


Sin(a)+^)T/2-] 
(a)+^)T/2  " J 
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-jajT/4n  ,T,  . 


T  T 

Because  of  the  relationship  between  and  (:*))#  it  is 


clear  that 


T  T 

lii  2  ( “  )  'I'  ) 


T  T 

i)2  (“w)  ’4)2  (w) 


so  that 


2 

'c ,  1 


=  o 


C  ,  2 


1  2 
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Thus  indeed,  a  =  c.  Observe  also  by  similiar  arguments,  th 


T  T 


I  q 


will  not  be  zero 


APPENDIX  B 


DETAILED  INVESTIGATION  OF  THE  BEHAVIOR  OF  THE  PRODUCTS  OF 

(-0))  AND  S'  (w) 

J  _ ^ 


We  have  defined 


F{S^  (t)p(t) } 


AT  c  _  ,  . 

-CO.)  + 

A  ^  J 


(B.l) 


Then 

F{S^(t)p(t)}  X  F{Sj;,(-)p(T)  }  =  S^(-w)Sj;(^) 

AT  AT 

=  -^ [ L  ( w  j  )  +  L  {-'.j+co^+w  j  )  ]  —  [L  ( )  +  L  ( +  )  j 

AT  ^ 

=  (“^)  [L(-5j-'jj  -co.)L(aj-(jj  -to,  )  +  L  ( -1.0+cj  +CJ  .  )  L  (  ) 

+  L  (  -  to  -  )  L  ( co+i,o^+tOj  )  +  L  ( -to'+!.o^+co  j  )  L  ( 'x'+.o^+ oj, )  ] 

(B.  2) 

Observe  that  for  reasonably  large  values  of  'jj^,  the  first  and 
the  last  term  in  this  expression  vanish,  and  we  are  left  with 
the  products 


The  orthogonality  condition  on  the  signals  required  that 
AwT  =  TT  or  Aa)T/2  =  ^7/2,  so  that  Eq.  B.5  becomes  (approximately) 


Eq.  B.6  is  zero  for  all  values  of  the  integer  j,  so  we  have 
that  for  AojT  =  the  product  S !  (-w)  Sj^  (w)  is  equal  to  zero  for 


P_^  ((ij)  =  /  Cos  co.te  dt  = 

Cl  0  ^  2 


^  3  T/2  Sin  (ujj^-u)  T/2 
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Furthermore 
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i  =  0,1 


The  third  term  in  Eq.  C,4  can  be  assumed  for  all  practical 
purposes  to  be  zero.  In  essence,  we  require  that  >>  2-/T 
i  =  0,1  for  the  approximation  to  be  correct. 

Consider  now  the  case  where 


S  i'il)  —  K[:('0  — LO-)  +  ‘‘{0  +  '..  •)] 
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then  K  =  ttP^  so  that  we  use 


S  (cj)  =  fr  p  [6(aj-(i).)  +  6  (oj+cj  .  )  ] 
c  c  3  D 


(C.  7) 


From  Eq.  C.2,  we  now  have 


2  2 
a  +  a  , 
c ,  o  c ,  1 


(C.  8) 


A  ssuming  that  co^  will  always  be  in  the  vicinity  of  and 
we  can  state  that 
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(uj,-<^.)T/2 


=  !Pel<--i> 


(C.ll) 


In  order  to  maximize  the  quantity  in  brackets  as  a  function  of 
uj j ,  we  need  to  take  derivatives  of  the  expression  and  set  it 
equal  to  zero.  The  result  of  this  operation  leads  to  a 
maximum  at  values  of  =  oj^  or  ojj  =  Therefore,  the 

maximum  value  becomes 


1  2  2  \ 

/  Sin  T/2 

\‘^c,o  *^0,1  /  max 

4  L 

)  +  1_ 

(C.13) 
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